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Abstract 

The presence of excess emission at 3.6-8.0 /iin was investigated in a sample of 27 
binary systems located in two nearby star-forming regions, Taurus and Ophiuchus, 
by using Spitzer /Infrared Array Camera (IRAC) archival data. Angular (Projected) 
separations for the binaries are greater than 2"(~280 AU), which allowed us to per¬ 
form spatially resolved photometry of individual primary and secondary sources. The 
measured occurrence of infrared excess suggests that binarity plays a role in the evolu¬ 
tion of circumstellar disks, even at such wide binary separations. Most of the binaries 
have excess emission from both the circumprimary and circumsecondary disks, or 
show photospheric levels for both components at all four wavelengths of IRAC. On 
the other hand, four systems (17lg 1 %, designated by “mixed” systems) exhibit excess 
emission from a single binary component. This ratio is significantly smaller than that 
predicted by the random pairing of single stars, suggesting that circumprimary and 
circumsecondary disks are synchronously dispersed. In addition, the excess frequen¬ 
cies (EFs) of primary and secondary sources with a projected distance of a p ~ 280- 
450 AU are 100t? 7 % and 9lt® 8 %, respectively, and significantly higher than that of 
single stars (70 ±5%). We made a simple model describing the EF distribution as 
a function of the disk outer radius, i? out . Comparisons with observations using the 
Kolmogorov-Smirnov test show that the observational data are consistent with the 
model when the EF ~ 1 region is found at R out ~ 30-100 AU. This disk radius is 
smaller than that typically estimated for single stars. The high EF of circumstellar 
disks with these radii may indicate a prolonged lifetime of dust in binary systems 
possibly because smaller disks counteract mass loss by photoevaporation. 

Key words: stars: binaries: general — stars: circumstellar matter — infrared: 
stars — stars: pre-main sequence 
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1. Introduction 


The ubiquity of binary stars is widely known. More than 60% of main-sequence solar- 
type stars (e.g., Abt & Levy 1976; Duquennoy & Mayor 1991) and ~30-40% of M-type stars 
have been reported as multiple systems (e.g., Fischer & Marcy 1992; Delfosse et al. 2004). In 
particular, the significant progress of the multiplicity survey of solar-type stars was achieved by 
Duquennoy & Mayor (1991) based on the radial velocity measurements. Raghavan et al. (2010) 
updated and extended their results using an increased sample with long-baseline interferometry 
and speckle interferometry. 

Numerous young binaries have also been detected in nearby star-forming regions (e.g., 
Ghez et al. 1993; Simon et al. 1993). Observational study of young binaries is particularly im¬ 
portant for understanding how most stars acquire their planetary systems. There are notable 
observations on circumstellar disks in binaries, although their evolutionary process is less clear 
than in single stars because sufficient data have not been accumulated. Circumprimary and cir- 
cumsecondary disks around SR 24 in Ophiuchus have been directly imaged in the near-infrared 
with the high spatial resolution, enabling a detailed comparison with theoretical predictions of 
disk formation and evolution in binary systems (Artymowicz & Lubow 1994; Mayarna et al. 
2010). On the other hand, excess emission in infrared and millimeter regions and line emis¬ 
sion as a diagnostic of accretion have been statistically analyzed to investigate the evolution 
of circumstellar disks in binaries (e.g., Harris et al. 2012). Recent ALMA observations yielded 
several new detections of low-mass circumsecondary disks thanks to the high sensitivity, and 
the increasing sample size helps to determine the disk population for each binary component 
(Akeson & Jensen 2014). 

In infrared and optical, consistent results have emerged from the statistical studies such 
as on the difference in disk frequency between primary and secondary stars, and the depen¬ 
dence of the disk lifetime on binary separation. White & Ghez (2001) investigated the nature 
of circumprimary and circumsecondary disks from their Ha and K — L colors. Because the 
lifetimes and accretion rates of circumprimary disks exceeded those of the secondary disks, 
they concluded that, in general, most of the material in a circumprimary disk is preferentially 
supplied by the circumbinary disk. McCabe et al. (2006) also studied the K — L, K — N colors, 
and Ha for binary systems. Their study showed that ~10% of the stars with the gas accretion 
and the excess in K — N color were photospheric in K — L color, indicating that their disks are 
in the middle of the inside-out evolution. These disks are harbored by the secondary stars with 
primary stars surrounded by disks, while only one of these is harbored by a primary star with a 
diskless secondary star; i.e., circumsecondary disks tend to disperse faster than circumprimary 
disks. Similarly, Monin et al. (2007) used Ha emission, K — L, and K — N to investigate circum¬ 
stellar disks in binaries. They found that there were fewer systems showing the gas accretion 
only onto the secondary stars than systems only accreting onto the primary stars. Therefore, 
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they deduced that the lifespan of circumprimary disks exceeds that of circumsecondary disks. 

Past studies have also revealed how the disk existence depends on binary separation 
(Prato & Simon 1997; Duchene et al. 1999; White & Ghez 2001; Hartigan & Kenyon 2003; 
Daemgen et ah 2012). Cieza et al. (2009) showed that the lifetime of circumstellar disks in 
binary systems was significantly shorter at separations of 40 AU and below compared to disks 
of single stars. They investigated infrared excess from binaries in four nearby star-forming 
regions (~l-3 Myr) using the [3.6] — [8.0] color obtained with Spitzer /IRAC without resolving 
the individual stellar components. Accelerated dispersal of circumstellar disks in close binary 
systems is also suspected among the members of p Chamaeleon, the ages of which are ~8 Myr 
(Bouwman et al. 2006). Similar results have been recently reported by Kraus et al. (2012), who 
studied excess emission in Taurus binaries over a wide range of wavelengths, from near-infrared 
to millimeter. From these data, they calculated the disk frequency, also without resolving 
individual stellar components. While ~67% of the close binaries (<40 AU) disperse their disks 
within ~1 Myr, ~80-90% of single stars and wide binary systems (>400 AU) retain their disks 
during 2.— 2 Myr. 

Some properties of the circumstellar disks, such as the disk frequency (e.g., Bouwman 
et al. 2006; Cieza et al. 2009; Kraus et al. 2012), gas accretion, and dust emission (e.g., White 
& Ghez 2001; Daemgen et al. 2012), have been reported to be indistinguishable between wide 
binaries and single stars. While many studies indicate rapid disk dispersal in close binary 
systems, the dependence of the disk dissipation time scale on wide binary separation (>100 AU) 
has been little investigated. Especially, sample size was insufficient, and binary components 
were not spatially resolved at mid-infrared wavelengths where the presence of disks can be more 
reliably estimated than in near-infrared. 

In this paper, we discuss infrared emission from each stellar component of binary sys¬ 
tems. The analysis is based on spatially resolved photometry at four wavelengths, 3.6, 4.5, 5.8, 
and 8.0 pm, obtained from Spitzer/TRAC images. The IRAC wavelengths correspond to a dust 
temperature of ~360-800 K and a radial region within a few AU for M-type stars. Thus, these 
wavelengths can probe the inner planet-forming radii in dusty disks. In addition, Spitzer/TRAC 
bands located between L and N that are difficult to access with ground-based telescopes; there¬ 
fore, new information such as the onset of infrared excess emission can be obtained. Moreover, 
the unprecedented sensitivity of Spitzer enables the detection of lower-mass secondaries. We 
collected binary systems with the projected separation a p of 2. ,/ 0-17."0, corresponding to ~280- 
2400 AU at 140 pc. Because the images of the primary and secondary stars overlap at projected 
separations near 2", we analyzed each star by stellar profile fitting for the photometry. 

The rest of the paper is organized as follows. Section 2 describes our sample and the 
Spitzer /IRAC data. Section 3 explains the profile fitting photometry . Section 4 presents the 
measured frequency of excess emission in binary systems and the difference between binary 
systems and single stars. The results are discussed in Sections 5 and summarized in 6. 
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2. Sample and archival data 

The full-width at half-maximum (FWHM) of a point source is 2" at 8.0 /mi for 
Spitzer /IRAC (Fazio et al. 2004). Given the array pixel scale of l."2, we selected binary systems 
with projected separations greater than 2"0. Thirty-three binaries of separation of H'95-1770 
in the star-forming regions of Taurus (distance 140 pc; Torres et al. 2007, age 1-2 Myr; Kraus 
& Hillcnbrand 2009a) and Ophiuchus (distance 140 pc; Loinard et al. 2008, age 1 Myr; Allers 
et al. 2006; h uh man & Rieke 1999; Prato et al. 2003) were extracted from the literature (Kraus 
et al. 2011; Kraus & Hillenbrand 2009b; Duchene et al. 2007; Ratzka et al. 2005; Haisch et al. 
2004; Haisch et al. 2002; White & Ghez 2001; Kohler & Leinert 1998; Reipurth & Zinnecker 
1993; Leinert et al. 1993) (see also Table 1 for references on each object). 

The companionship and membership in Taurus or Ophiuchus in our sample were ex¬ 
amined using the literature because the sample may include suspicious secondary stars not 
physically associated to the primary, or older systems than 1-2 Myr. In Taurus, coevality 
(dAlog(f)) between primary and secondary stars has been indicated to be <0.4 dex which is 
significantly smaller than 0.58 dex estimated from the random paring of single stars (Kraus & 
Hillenbrand 2009a). In our initial 33 binaries, the systems whose ages are known and coevality 
is larger than 0.4 dex are 2MASS J04554757+3028077, HP Tau/G2, and DK Tau. The age dif¬ 
ference of dAlog(f) ~ 0.8 was found for DK Tau between the primary and the secondary stars, 
but a possible systematic error was implied in this age estimation (Kraus & Hillenbrand 2009a). 
In addition, in order to explain the estimated size for the circumprimary disk, the presence of 
a companion is favored just by the distance of the secondary candidate. Therefore, we retained 
DK Tau in our binary sample while excluded 2MASS J04554757+3028077 and HP Tau/G2. In 
addition, RXJ 0437.4+1851, NTTS 040142+2150, and V1117 Tau systems are speculated to 
be older compared to the typical ages of Taurus based on the Li test, or to be unassociated to 
the Taurus region judging from the proper motions (Sestito et al. 2008; Wichmann et al. 2000). 
Since contamination of these sources can affect the discussion of disk frequency at the age of 
1-2 Myr, we eliminated these sources from our sample. For Ophiuchus, VSSG 23 was removed 
from our sample due to their non-coevality (dAlog(f) ~ 1; Prato et al. 2003). LInfortunately, 
sufficient data on coevality and proper motions were not available for many of the Ophiuchus 
sources, and these systems remain to be included in the sample as binary candidates. As a 
result, six systems were eliminated from the initial 33 binaries, and 27 systems were analyzed 
in this paper. 

It should be noted that our sample includes triple and quadruple systems consisting of 
the close binaries with separations of ~0701-0f'l and the other widely separated components: 
2E 1628.2-2423 (Mathieu et al. 1989), V1001 Tau (White & Ghez 2001), RXJ 0435.9+2352 
(Kohler & Leinert 1998), L1689 SN02 (Ratzka et al. 2005), and 2MASS J04251767+2617504 
(Duchene 1999). Because the close components should exhibit rapid disk clearing (Monin et 
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al. 2007; Cieza et al. 2009; Kraus et al. 2012), whether this parameter affects our discussion 
or not is explored in later sections. Binarity in these systems, especially those that are widely 
separated, has yet to be confirmed by multi-epoch astrometry; therefore, many of these systems 
are candidates of binaries. Table 1 summarizes the projected separations and spectral-types 
taken from literature for our sampled systems. 

Since the closely spaced systems cannot be analyzed by aperture photometry, they were 
divided into individual components by profile fitting to the Spitzer /IRAC archival images, as 
described in the following section. Fortunately, many of the binaries are located in regions 
surveyed in the “From Molecular Cores to Planet-Forming Disks” (c2d) legacy project (Evans 
et al. 2003). The IRAC data were acquired from the NASA/IPAC Infrared Science Archives 1 . 
Following basic calibration, cosmic rays were removed by mosaicking the post-BCD data (IRAC 
Instrument Handbook 2 ), and photometry was performed on the mosaicked post-BCD images. 
The pixel scale of the post-BCD images was 0."6, two times as small as the native pixel size. The 
data were obtained over two effective integration times: long (10.4 sec) and short (0.4 sec). 
Photometry was generally carried out using the short-exposure frames, because these frames 
yield better fits to the stellar profile. The exceptions were the faint (<9 mag) binaries for which 
the long-exposure images were used. 

3. PSF-fitting photometry 

The images of the primary and secondary stars overlapped at binary projected separa¬ 
tions of ~2 ,/ ~4 // . Therefore, we resolved the primary and secondary stars by y 2 fitting using a 
stellar profile model, which determined six variables: the positions of the primary and secondary 
star (x and y for each star) and the scaling factors of the profile models. The six simultaneous 
equations were solved by the Newton-Raphson method. The initial position of the primary 
star was set at the location of the brightest pixel and was adjusted to minimize y 2 . The lo¬ 
cation of the secondary star was assumed to match its relative position to the primary star, 
reported in the literature. The location of the secondary star was then adjusted to minimize 
y 2 . The initial scaling factor for the profile model was determined to match the pixel intensity 
to its observed intensity at the initial stellar position. To create the profile model, the images 
of Eve single stars near the binary system were averaged using the PSF task in IRAF 3 . The 
uncertainty on the observed PSF used for the y 2 fitting was estimated as follows. First, each of 
the five single stars was fitted with the same fitting tool and profile model, using the standard 
deviation o stored in the 2D uncertainty image in the post-BCD dataset (*_unc.fits) as a 

1 http://irsa.ipac.caltech.edu/data/SPITZER/docs/spitzerdataarchives 

2 http://irsa.ipac.caltech.edu/data/SPITZER/docs/irac/iracinstrumenthandbook/ 

3 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association 
of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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weight. Here, a contained dark noise, flat-fielding error, read-out noise, and Poisson noise. 
At each pixel, the difference between the observed value and the fitted model profile reached 
~ 10-20%, exceeding those predicted from the uncertainty image by a few orders of magnitude. 
Thus, to further consider the possible uncertainty associated with the observed PSF, the mean 
of the residuals between the observed profile and the model was calculated from each of the five 
stars, then applied to the y 2 fitting as a new a. Averaged stellar profile model, together with 
its a, yielded plausible results (y 2 /d.o.f. ~ 1). We also confirmed the quality of the fitting in 
the residual image and that the location of each star was consistent in the different wavebands, 
even for closely spaced systems (~2 // -3 ,/ ). The y 2 of all targets are summarized in Table 1. 
The widely separated systems (>4") were fitted with the same tool to prevent any systematic 
error introduced by the fitting tool. 

The averaged stellar profile could not be obtained for one system (TYC 1289-513-1), 
because no nearby single stars were found in the same field. This system was fitted by the 
template point response function (PRF) provided by the IRAC team. Since the FWHM of the 
template function was always smaller than that of the observed PSF, it was convoluted with a 
Gaussian function. The values store in the post-BCD uncertainty image was used as a for the 
y 2 fitting. The resulting y 2 /d.o.f.(> 1000) were unacceptably large, which should be plausible 
because the uncertainty associated with the observed PSF was ignored. We verified the effect 
of variation in the stellar profile on the photometric results using systems that permitted an 
averaged stellar profile. The magnitudes measured between the Gaussian-convolved template 
and averaged stellar profile fitted to the data differed by <0.1 mag in all the IRAC wavebands. 
Therefore, we concluded that the cr in the uncertainty files were excessively small and responsible 
for the huge y 2 , but the Gaussian-convolved template was reasonably fitted to the data. 

The photometric error can be estimated using the revised a obtained from the empirical 
PSF constructed from the stellar profiles of five nearby single stars. The best-fitting scaling 
factor of the profile was varied until the reduced-y 2 changed by 1.17 (d.o.f. = 6), and the flux 
difference created by this process was considered to be an error in the measured magnitude. 
The uncertainty estimates of several binaries were consistent within 10% in all four wavebands. 
Comparing the measured magnitudes with those obtained by Gaussian-convolved template 
fitting, we estimated a conservative photometric error of 0.14 mag in all the IRAC bands. 

Photometry in some IRAC bands was precluded in three systems, because the primary 
and secondary stars were obscured by a large brightness contrast (at 3.6 /mi for DoAr 24E and 
at 8.0 and 5.8 //in for HN Tau B) and by the overlap between the secondary position and the 
latent image caused by an exceptionally bright primary star (at 8.0 /irn for 2E 1628.2-2423). 
The measured magnitudes of all targets are summarized in Table 1. 
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Table 1 . Characteristics of Sampled Binary Systems. 


Systems 

Separation 

SpT. 

[3.6] 

[4.5] 

[5.8] 

[8.0] 

X 2 /d.o.f. (chl-4) 


(”) 


(mag) 

(mag) 

(mag) 

(mag) 


Taurus 

JH 223 

2.07 (3) 

M2 (11) 

9.10 

8.88 

8.39 

7.88 

1.0, 0.4, 0.7, 1.0 



M6.5 (11) 

11.26 

11.06 

10.63 

9.85 


CoKu Tau 3 

2.1 (14) 

Ml (11) 

7.51 

6.98 

6.51 

5.67 

0.8, 1.0, 1.5, 1.6 



M5 (11) 

9.17 

8.99 

8.50 

7.99 


DK Tau 

2.3 (2) 

K9 (9) 

6.21 

5.77 

5.50 

4.87 

1.0, 1.2, 1.2, 1.3 



Ml (9) 

7.84 

7.47 

7.26 

6.72 


CIDA 9 

2.34 (3) 

MO (11) 

9.33 

8.74 

8.09 

7.08 

1.6, 1.0, 0.9, 1.3 



M2.5 (11) 

11.49 

11.01 

10.70 

9.82 


HK Tau 

2.34 (3) 

Ml (1) 

7.81 

7.49 

7.14 

6.62 

0.7, 0.8, 0.4, 0.8 



Ml (1) 

11.47 

10.91 

10.64 

9.69 


IT Tau 

2.39 (2) 

K3 (22) 

7.54 

7.16 

6.82 

6.32 

1.8, 1.9, 0.8, 0.7 



M4 (22) 

9.05 

8.68 

8.11 

7.51 


V1001 Taut 

2.62 (2) 

K8 (9) 

6.69 

6.27 

5.88 

5.33 

0.8, 1.8, 1.1, 0.9 



MO (9) 

7.90 

7.38 

7.11 

6.06 


HN Tau 

3.10 (6) 

K5 (4) 

7.00 

6.34 

5.65 

4.76 

1.3, 3.0,—,— 



M4.5 (22) 

11.11 

10.51 

— 

— 


V710 Tau* 

3.17 (2) 

Ml (21) 

7.99 

7.80 

7.19 

6.57 

1.6, 4.1, 1.4, 3.1 



M3 (21) 

8.39 

8.45 

8.26 

8.29 


2MASS J04251767+2617504t 

3.4 (2) 

K7 (13) 

8.44 

8.38 

8.34 

8.31 

2.9, 1.2, 1.8, 1.3 



M3 (13) 

9.10 

9.05 

8.98 

9.00 


BBM92 14 

4.06 (12) 

M2 (10) 

6.55 

5.84 

5.32 

4.57 

1.2, 1.2, 1.6, 2.2 



M2 (10) 

7.13 

6.53 

5.97 

5.14 


TYC 1289-513-1 

6.87 (13) 

— 

8.68 

8.67 

8.57 

8.47 

— 



— 

10.69 

10.69 

10.58 

10.47 


RXJ0435.9+2352t 

11.3 (13) 

— 

8.79 

8.78 

8.76 

8.63 

9.2, 5.3, 2.3, 0.8 



— 

10.82 

10.98 

10.84 

10.69 


GI Tau 

13.14 (6) 

K6 (21) 

6.43 

6.14 

5.90 

4.90 

0.8, 1.0, 1.8, 1.0 



M0 (21) 

6.83 

6.38 

5.81 

4.82 


FZ Tau 

17.17 (6) 

M0 (4) 

6.28 

5.75 

5.28 

4.62 

1.5, 0.7, 1.6, 1.7 



K7 (21) 

7.19 

6.86 

6.53 

6.07 


Ophiuchus 

DoAr 24 E 

2.07 (15) 

G6 (7) 

— 

5.03 

4.62 

3.88 

—,1.8, 2.1, 1.4 



K4.5 (7) 

— 

6.25 

5.77 

4.93 
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Table 1. (Continued.) 


Systems 

Separation 

(”) 

SpT. 

[3.6] 

(mag) 

[4.5] 

(mag) 

[5.8] 

(mag) 

[8.0] 

(mag) 

X 2 /d.o.f. (chl-4) 

L1689 SN02 t 

3.01 (15) 

M2 (10) 

6.34 

5.50 

4.76 

3.96 

— 



M2 (10) 

8.04 

7.72 

7.37 

6.76 


WSB71a 

3.56 (15) 

K2 (5) 

6.62 

6.08 

5.46 

4.75 

1.7, 0.7, 2.4, 1.3 



M6 (5) 

9.49 

9.15 

8.53 

7.80 


ROXRF 36 

3.59 (15) 

— 

9.26 

9.16 

8.13 

9.01 

1.2, 0.8,—, 0.2 



— 

10.09 

9.96 

9.80 

9.77 


WL 18 

3.62 (15) 

— 

8.67 

8.00 

7.84 

7.35 

2.0, 2.7, 1.5, 0.3 



— 

11.06 

10.68 

10.34 

10.22 


2MASS J16264848-2428389 

4.15 (16) 

— 

9.08 

8.21 

7.56 

6.74 

1.9, 3.3, 1.4, 1.2 



— 

10.38 

9.71 

8.94 

8.25 


2E 1628.2-2423^* 

4.3 (14) 

GO (20) 

6.37 

5.96 

5.30 

— 

2.0, 1.0, 2.7,— 



K5 (20) 

7.68 

7.59 

7.50 

— 


2MASS J16262097-2408518* 

5.21 (15) 

M3 (4) 

9.14 

9.13 

9.06 

8.98 

2.2, 1.0, 0.4, 0.2 



M7 (4) 

10.89 

10.53 

10.10 

9.37 


2MASS J16271757-2428562 

8.55 (16) 

F7 (18) 

7.82 

6.95 

6.47 

6.49 

5.2, 7.8, 6.0, 2.0 



M4 (19) 

8.36 

7.35 

6.51 

5.95 


2MASS J16263682-2415518 

9.08 (12) 

M0 (3) 

8.17 

7.85 

6.29 

5.27 

4.1, 2.3, 2.7, 1.5 



M5 (7) 

10.25 

10.23 

9.74 

9.14 


UCAC2 21797671* 

10.3 (14) 

M2 (7) 

8.83 

8.48 

8.23 

7.72 

0.8, 1.9, 1.7, 0.8 



M4.5 (7) 

10.07 

10.04 

9.96 

9.87 


2MASS J16262404-2424480 

10.47 (17) 

— 

6.43 

4.83 

4.18 

3.44 

1.6, 5.0, 2.9, 3.8 



— 

7.44 

6.55 

5.78 

5.04 



NOTE.—For individual binary systems, spectral type (Column 3) and measured IRAC magnitudes (Columns 4-7) of the primary 
and secondary members are shown in the first and second row, respectively. The second column lists the projected separations of 
the systems. The eighth column lists the x 2 /d.o.f- given by the profile fitting at each IRAC wavelength. References are given in 
parentheses in Columns 2-3 (see below). The typical photometric error in the PSF fitting photometry is 0.14 mag. 

* Denotes a candidate of a triple or quadruple system. 

Denotes a mixed system. 

(1) White Sz Hillenbrand 2004; (2) White Sz Ghez 2001; (3) Kraus et al. 2011; (4) Correia et al. 2006; (5) Prato et al. 2003; (6) 
Kraus Sz Hillenbrand 2009b;(7) Wilking et al. 2005; (8) Kraus Sz Hillenbrand 2009a; (9) Monin et al. 1998; (10) Connelley Sz Greene 
2010; (11) Kraus et al. 2007; (12) Duchene et al. 2007; (13) Kohler Sz Leinert 1998; (14) Reipurth Sz Zinnecker 1993; (15) Ratzka 
et al. 2005; (16) Haisch et al. 2004; (17) Haisch et al. 2002; (18) Greene Sz Meyer 1995; (19) Luhman Sz Rieke 1999; (20) Walter et 
al. 1994; (21) Kenyon Sz Hartmann 1995; (22) Duchene et al. 1999. 







4. Results 


4-1. Excess emission in IRAC bands 

We first examined whether the emission of the individual components in the binary 
systems exceeds the photospheric level at 8.0 /an. Figure 1 shows a histogram of the color 

[3.6] — [8.0] distribution of the sources for which the photometry at 3.6 and 8.0 /irn is available. 
The clear gap in the figure, occurring at [3.6] — [8.0] ~ 0.6, is considered to divide the sources 
with and without excess at 8.0 /irn. The excess frequency (defined as (excess sources)/(excess 
sources + non-excess sources), and hereafter denoted EF) of the primary and secondary stars 
were 79i7%(19/24) and 73jlf%(16/22), respectively. 



Fig. 1. Histogram of the color of [3.6] — [8.0] shown by primary ( filled) and secondary ( hollow ) sources. 

Class I sources (see Section 4.2) are not included in this plot. 

Figure 2 shows color-color diagrams of [3.6] — [4.5] (left) and [3.6] — [5.8] (right) against 

[3.6] — [8.0]. The objects with [3.6] — [8.0] < 0.6 are gathering near the origin of the diagrams, 
indicating that they are also photospheric at 4.5 and 5.8 /irn. In contrast, the sources with 

[3.6] — [4.5] > 0.2 and [3.6] — [5.8] > 0.4, located in the correlation extending to the upper right 
part of the diagrams, tend to have excess emission at 8.0 /irn ([3.6] — [8.0] > 0.6) in addition to 
4.5 and 5.8 /mi. However, it should be noted that identification of the excess sources at 4.5 /irn 
is difficult because no clear color break occurs. Although we can see a clear gap at [3.6] — [5.8] 
color around 0.4, some excess sources overlap with the non-excess sources when considering the 
typical uncertainty of the [3.6] — [5.8] colors. 

The 8.0 /iin excess was detected relative to 3.6 /mi, rather than absolutely, based on 
the color [3.6] — [8.0]. In addition, the excess at 4.5 and 5.8 /irn is sometimes difficult of 
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Fig. 2. Color-color diagrams of [3.6] — [4.5] (left) and [3.6] — [5.8] (right) against [3.6] — [8.0]. Filled circles 
show primary stars and hollow circles indicate secondary stars. The extinction vector in the IRAC bands 
is based on the simple-fit formula reported by Indebetouw et al. (2005). Class I sources (see Section 4.2) 
are not included in this plot. 

identification by the color-color diagrams. Therefore, we additionally checked the presence of 
excess emission using the spectral energy distributions (SEDs). Fitting a stellar photosphere 
requires photometry in one or more optical bands in addition to J band. The photometric 
values in the optical bands (B, V, R , I ) were taken from the USNO-B1 Catalog (Monet 2003). 
For the sources for which the resolved photometric values were unavailable from the catalog, we 
used the resolved data in optical bands from additional references (Zacharias et al. 2005; Torres 
et al. 2006; Grankin et al. 2007; Herbig & Bell 1988). The J-band magnitudes were taken 
from the 2MASS All-Sky Catalog of Point Sources (Cutri et al. 2003). If these magnitudes 
were unavailable from the catalog, the resolved data in the J band were extracted again from 
additional references (Kraus et al. 2007; Connclley & Greene 2010; Moneti & Zinnecker 1991; 
Chavarria et al. 2000). The spectral types were available for many of our sample sources in 
previous studies (see Table 1). Free parameters for the fitting are Ay and the scaling factor 
of the photosphere, corresponding to the radius of the star and distance to the star-forming 
regions. To model the stellar photosphere at high and low effective temperatures, we adopted 
the models of Krucz (Kurucz 1993) and AMES-Cond (Allard et al. 2001), respectively. For 
sources with known Ay but unknown optical-band magnitudes, the quoted Ay values were used 
to correct for the extinction in the SEDs (Furlan et al. 2011; Chavarria et al. 2000; Wahhaj et 
al. 2010; Evans et al. 2009; Currie & Sicilia-Aguilar 2011; McClure et al. 2010; Bontemps et al. 
2001). The data required for the SED fitting were unavailable for 12 of the sources, and we do 
not discuss them below. 

The excess at 8.0 /.irn was estimated based on the SEDs for 37 sources with the available 
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photometry at this wavelength. All the sources which appear at [3.6] — [8.0] > 0.6 show the 
excess at 8.0 /an in their SEDs, while none of them at [3.6] — [8.0] < 0.6 exhibit an excess 
emission. We also confirmed that 34 out of the 37 sources (92±g%) showed or did not show 
excess at both 3.6 and 8.0 /im. Six sources (ROXRF 36 A and B, TYC 1289-513-1 A and 
B, RXJ0435.9+2352 A and B) were excluded in the above calculation, because they lacked 
spectral-type data. However, these sources were assumed to be non-excess sources in all the 
IRAC bands, since they were classified as Class 111 based on A" — [8.0] data (see next subsection). 
If the six sources are included, 40 out of the 43 selected sources (93±g%) either show or do not 
show excess at both 3.6 and 8.0 fin r. 

The occurrence of excess emission and the relation with the 8.0 /irn excess were also 
examined at 4.5 /im and 5.8 fim in the SEDs. Thirty nine out of 54 sources were used since 
the photometric data were obtained for them at these wavelengths and 8.0 /im. Note that the 
sample for 5.8 /im is different from that of 4.5 /im, although the sample size is the same. At 
4.5 /im, 37 out of the 39 sources (95ig%) showed excess or non-excess at both wavelengths 
(4.5 /im and 8.0 /im), and at 5.8 /im, 38/39 sources (98±4%) present the same characteristics. 
These results indicate that the presence or non-presence of excess emission at 8.0 /an behaves 
similarly at 3.6, 4.5, and 5.8 /im. The exceptions are three sources that show excess at 8.0 /an 
but not at the shorter wavebands. 2MASS J16263682-2415518 B, V 710 Tau A, and JH 223 A 
begin to exhibit significant excess from ~8.0, 5.8, and 4.5 /an, respectively. These three sources 
are further discussed in the appendix. 

4-2. Our sample of disk-bearing stars 

Following the results in the previous subsection, we evaluate the existence of a circum- 
stellar disk through the color of [3.6] — [8.0] since it shows the most distinct boundary between 
excess and non-excess stars. However, the magnitudes at 8.0 /an were unavailable for two 
sources (2E 1628.2-2423 and HN Tau B); thus we attempted to speculate whether the sources 
lacking in the photometric data at 8.0 /an had excess at this wavelength (i.e., [3.6] — [8.0] > 0.6). 
2E 1628.2-2423 A shows significant excess at 3.6-5.8 /an, as confirmed from the SED, thus it 
is plausible to categorize it as an 8.0 /an excess source. On the contrary, 2E 1628.2-2423 B 
exhibited no excess at 3.6-5.8 /an. Considering the rarity of sources that first show excess 
at 8 /an (~6%), this object was presumed as a non-excess source at 8.0 /an. The value of 
[3.6] — [4.5] = 0.41 for HN Tau B (M4.5; Kraus & Hillcnbrand 2009a) is significantly higher 
than the photospheric color of an M5-type star (—0.08; Luhman et al. 2010). Therefore, this 
source most likely has 8.0 /tm excess. Note that the recent ALMA observations failed to detect 
the secondary source in dust continuum at 850 /an and 1.3 mm (Akeson & Jensen 2014), but 
the secondary indicated the Ha emission and was classified as a classical T Tauri star (e.g., 
Duchene 1999). 

Photometry at 3.6 /an was not obtained for both components in DoAr 24 E (see 
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Section 3), and hence the [3.6] — [8.0] colors were unavailable. However, they generated the 
significant excess at 4.5-8.0 /an, as confirmed by the SEDs. Therefore, they can be catego¬ 
rized as excess sources. In summary, in the statistical discussion in later sections, we regard 
2E 1628.2-2423 A, HN Tau B, and DoAr 24 E as excess sources and 2E 1628.2-2423 B as a 
non-excess star. 

Our sample includes Class I sources that might still be surrounded by envelopes. To 
simplify the discussion and to calculate the excess frequency for Class II and III sources only, we 
excluded the Class I sources as follows. Classes were identified from the A" — [8.0] color, adopting 
the criterion of Lada (1987). That is, sources satisfying K — [8.0] > 3.9 were regarded as Class I, 
1.14 < K — [8.0] < 3.9 as Class II, and —0.26 < K — [8.0] < 1.14 as Class III. The JH A'-bands 
magnitudes were taken from the 2MASS Catalog or additional references (Kraus et al. 2007; 
Connelley & Greene 2010; Moneti & Zinnecker 1991; Chavarria et al. 2000). For the sources 
lacking in the SED fitting (Section 4.1), the Ay values for the color correlation were taken from 
previous studies (Furlan et al. 2011; Kraus & Hillenbrand 2009a; Wahhaj et al. 2010; Chavarria 
et al. 2000; Evans et al. 2009; McClure et al. 2010; Currie & Sicilia-Aguilar 2011; Bontemps 
et al. 2001). As a result, 35 out of the 54 sources with both K — [8.0] and Ay data available 
were divided into their respective classes, and two sources (2MASS J16271757-2428562 B and 
CIDA 9 A) were identified as Class I. The remaining 19 sources lacking K — [8.0] or Ay data 
were classified based on the J — H versus H — K color-color diagram. In this classification, we 
employed the JHK colors of unreddened main sequence of Koornneef (1983) (Lada & Adams 
1992) and the locus of T Tauri stars derived by Meyer et al. (1997). In addition, we assumed the 
interstellar reddening raw estimated by Koornneef (1983), Bessell & Brett (1988), and Martin & 
Whittet (1990). As a result, one source (2MASS J16262404-2424480 B) was classified as Class I. 
The 51 sources or 24 systems not classified as Class I were used for the analyses described in 
Section 5. 

5. Effect of binarity for IRAC excess and differences from single stars 

In this section, we show whether primary and secondary stars exhibit the same behavior 
regarding the presence of IRAC excess. In addition, we discuss the dependence of EF on binary 
separation and the difference of the EF between the binary systems and single stars. 

5.1. Mixed systems 

In Figure 3, the [3.6] — [8.0] color of the secondary star is plotted against that of the 
primary star. As described in the previous section, the clear color break at ~0.6 divides the 
sources showing excess from those not showing excess. Twenty systems out of the 24 binaries 
(831?!%) show excess or non-excess associated with both components. The remainder (17lg 1 %) 
show excess emission from only one of the binary components (“mixed” binaries). The number 
of systems showing excess solely from the primary and that solely from secondary source are 
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Fig. 3. [3.6] — [8.0] color of secondary versus primary stars. The extinction vector in the IR.AC bands is 
based on the simple-fit formula reported by Indebetouw et al. (2005). The typical error in the measured 
color is 0.20 mag. The sources without measurements of [3.6] — [8.0] for both binary components are not 
included in the figure. 

three (13/g%) and one ( 4 / 4 %), respectively. To compare this mixed system ratio with the 
distribution of random sampling from a single star population, the EF of single stars was set to 
70% (see Section 5.3). The probability of obtaining the observed number of mixed systems is 
C|4(0.42) 4 (0.58) 20 ~ 6 x 10~ 3 . This significantly low probability implies that the disappearance 
of excess emission from circumprimary and circumsecondary disks is strongly correlated. On 
the other hand, a binomial statistical analysis revealed no preference for which component is 
exhibiting excess; i.e., there is no significant difference between the ratio of systems showing 
excess solely from the primary source (3/4; 75/ 2 g%) and from the secondary source (1/4; 
25 / 20 %)- The errors stated here represent the la confidence interval containing the central 
68% of the binomial distribution. 

Our sample includes five triple or quadruple systems with close (< 0".l ~ 15 AU) binaries 
corresponding to the primary or secondary component. Binaries with separations closer than 
~40 AU show a lower frequency of infrared excess compared with single or wide binary sources 
(e.g., Monin et al. 2007; Cieza et al. 2009; Kraus et al. 2012). This suggests that these triples and 
quadruples should have a higher fraction of mixed pairs. However, the EF does not significantly 
change when these five systems are excluded from the analysis. In this case, in 16 out of 19 
systems ( 84 / 43 %), both members either show or do not show excess emission. The number of 
mixed systems is 3 (16/g 3 % of the sample), and the conclusion is unaltered. 

The four mixed systems are V710 Tau, 2MASS J16262097-2408518, 2E 1628.2-2423, 
and UCAC2 217971. They can be in the evolutionary transition between disk-bearing and 
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disk-less binaries. Based on the observed ratio of mixed systems and the ages of Taurus and 
Ophiuchus which are ~10 6 yrs, the duration of the mixed phase can be speculated as <10 5 yrs. 
Alternatively, as mentioned above, it is known that the existence of a tight companion leads to 
a shorter timescale for disk clearing. Since very high angular-resolution observations have not 
been reported for the 4 binaries, we cannot exclude the possibility of these systems possessing 
undetected close (<10 AU) companions. 

The rarity of mixed systems contradicts the concept of random pairing suggested in 
previous studies. In their study of circumstellar disks in the Taurus region, White & Ghez 
(2001) reported that the disks at the separations wider than 210 AU are consistent with random 
pairing. This conclusion was inferred from the spatially resolved K — L color and Ha data. More 
recently, Daemgen et al. (2012) suggested that random pairing occurs in the systems separated 
by more than 200 AU in the Orion Nebula Cluster. In contrast, our analysis indicates the 
correlation of primary and secondary components in the separation range of ~280-2400 AU. 
Different sample size may account for the disagreement since our sample of 24 systems is 
approximately three times greater than those of the previous works. Another potential source 
of disagreement is the disk indicators. In the above-mentioned studies, gas accretion and dust 
emission were probed from the innermost regions. By contrast, the [3.6] — [8.0] color traces dust 
emission from the farther, ~1 AU region of the disk. 

5.2. Dependency of EF on separation 

The [3.6] — [8.0] color is plotted as a function of binary projected separation (a p ) in 
Figure 4. The horizontal axis is converted from arc seconds to AU, assuming that each star¬ 
forming region is located at 140 pc. In this figure, more closely separated binary companions 
(a p ~ 280-450 AU) tend to show excess emission. We verified the dependence of the EF on 
binary separation and distinguished the regions exhibiting high and low EFs employing Fisher’s 
exact test below. 

To identify the different EF regions, we split our sample into two groups based on 
whether the projected separation was larger or smaller than the separation a c and assessed 
whether EF differed between the groups. Figure 5 shows the results of Fisher’s test applied 
to the three sub-populations; primary, secondary, and both components. The left and right 
panels show Fisher’s p-value and EF, respectively, as a functions of a c . In the left panel, the 
p -value initially decreases at smaller a c , and is minimized at a c ~ 450 AU. This is because the 
sample size with a p < a c increases while maintaining the high EF, whereas that of binaries with 
a p > a c remains at ~75% as shown in the right panel. Fisher’s exact test clarified that the EFs 
of more closely and widely separated systems are different with a boundary at ~450 AU, at the 
significance of 1-2 <j for each binary component. Including both primary and secondary sources 
in the analysis, Fisher’s test gave the significance of 2-3 a. It is notable that the EF is very 
high (~100%) at a p < 450 AU for each primary and secondary subpopulation, and naturally 
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Fig. 4. [3.6] — [8.0] versus binary separation for primary (black circles) and secondary (unfilled circles ) 
stars. Separation was converted to physical distance (in AU) using the distance to each star-forming 
region (140 pc). Note that the separation is the projected separation. The sources of excess emission are 
clustered around 280-450 AU. The sources without measurements of [3.6] — [8.0] are not included in the 
figure. 

for the sample including both components (right panel of Figure 5). 

5.3. Comparison of binary systems and single stars 

To compare the EFs between the binary systems with a p < 450 AU and single stars, we 
first estimated the EF of single stars (EF S i ng i e ) in Taurus and Ophiuchus. In Taurus, the sample 
of 39 stars was adopted from Kraus et al. (2012) because little contamination of tight binaries is 
expected given the unprecedented resolution (~ 15-20 mas; 2-3 AU) and sensitivity (A K ~5-6 
mag at 40 mas;~ 7-15 Mj at ~6 AU) in their observations. The single stars in Ophiuchus were 
collected from Cieza et al. (2009). Unfortunately, the single nature of these stars is less clear 
compared to the Taurus stars and the sample can be contaminated by tight binaries. Since 
there is a possibility of EF lowered for single stars because of this contamination, we attempt to 
calibrate the EF for Ophiuchus stars based on the frequency of tight binaries found in Taurus, 
in the discussion later in this section. 

The EF can be biased by stellar age. In fact, the disk frequency is well-correlated with 
the typical age of a star-forming region (Mamajek 2009). Theoretically, the evolution of a 
pre-main sequence star at ~1 Myr with a spectral type of G-M is characterized by a decrease 
in luminosity at a nearly constant temperature (Baraffe et al. 1998; Yi et al. 2001). This means 
that the stellar luminosity can be substituted for the age; that is, the higher the luminosity, 
the younger the star and vice versa. Thus, we checked whether the age dependence of EF sing i e 
existed in our sample through the luminosities, which were estimated from the SEDs. Sufficient 
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Fig. 5. This figure shows the significance of EF differences between closer ( a p < a c ) and wider ( a p > a c ) 
binary systems. Left: Fisher’s p-value vs. a c plane. The differences in EF were tested by Fisher’s exact test 
between two groups of the binary samples, divided by whether the projected separation exceeds a c . The 
solid , dashed , and dotted line show the Fisher’s p-value for primary, secondary, and primary+secondary 
sources, respectively. The horizontal-dotted and dashed-dotted line indicate the 2 a and 3cr levels, re¬ 
spectively. The minimum p-value occurs at ^450 AU. At this radius, the significance of the EF difference 
between a p < a c and a p > a c is maximized (2er <p< 3cr for each component and >3<r for primary+secondary 
sources). Right: The EF of each group bifurcates in all three sub-populations. 

data for the SED fitting were available for 31 Taurus and 62 Ophiuchus sources. 

To examine the dependence of EF on luminosity, the Taurus sample was divided into 
two subpopulations of spectral types of G0-K7 and M0-M4, which are comprised of 12 and 19 
sources, respectively. Each subpopulation was then divided by luminosity into two groups of 
approximately the same size. The EF S i ng i e was determined from the [3.6] — [8.0] colors. The mag¬ 
nitudes at these wavelengths were obtained from IPAC/Gator 4 . When the data were unavail¬ 
able, the photometric results were extracted from the literature (Wahhaj et ah 2010; Luhman 
et al. 2010). A systematic difference between the above-cited magnitudes and those measured 
with our PSF-fitting tool do not affect the EF comparison between single and binary systems, 
because in both cases, the [3.6] — [8.0] distribution is clearly separated into two groups: with 
and without the color excess. Using the same criteria ([3.6] — [8.0] = 0.6), the EFs of four sub¬ 
populations were indistinguishable (EF~0.7-0.9), indicating that EF S i ng i e does not significantly 
depend on luminosity. Thus, we infer that the EF sing i e is independent of stellar age. When 
we added seven M0-M4 stars without sufficient observations for SED fitting, the 38 sources 
of Taurus gave an EF of 74% (28/38). The same analysis was applied to the Ophiuchus sam- 

4 http://irsa.ipac.caltech.edu/applications/Gator/ 
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Fig. 6. This figure shows the significance of EF differences between closer binary systems ( a p < a c ) and 
single stars without the calibration on possible tight companions described in the text. Left: Fisher’s 
p- value vs. a c plane for primary only, secondary only, and primary+secondary sources (solid, dashed , and 
dotted lines, respectively). The EF difference between the binary sample with a p < a c and single stars 
was assessed by Fisher’s exact test. The horizontal-dotted and dashed-dotted lines indicate the 2a and 3cr 
levels, respectively. As in Fig. 5, the p- value is minimized at ~450 AU, where the difference in EF between 
binary systems with a p < a c and single stars becomes most significant (>2er for each component and >3er 
for primary+secondary sources). Right: The EF of binary systems and single stars. The horizontal-dotted 
line indicates the EF s j ng i e (63%). 

pie, and a similar result was obtained about the EF dependence on the stellar luminosity. To 
equalize the comparison between the single and binary samples, we discarded one object of the 
spectral type later than M6 because our binary sample consists of <M6 members. The EF was 
calculated as 56% (34/61) for the Ophiuchus cample. The final sample of single stars consists 
of 38 in Taurus and 61 sources in the Ophiuchus regions. 

The difference in EF between single and binary systems was investigated first using this 
whole sample of single stars. Out of 99 sources, 62 show excess ([3.6] — [8.0] > 0.6); i.e., EF sing i e 
is 63% (62/99). This EF sing i e was then compared with the EF of binaries with separations of 
a p <a c , again using Fisher’s exact test. Figure 6 shows the results for the three subpopulations 
(primary, secondary, and both). The left and right panels present the Fisher’s p -value and EF, 
respectively, as functions of a c . The most significant difference appears at a c ~450 AU, because 
the EF of primary and secondary sources separated within ~450 AU is unity as described in 
Section 5.2. This EF of each component is significantly (>2er) higher than that of single stars of 
63%. For the entire binary sample including both primary and secondary sources, Fisher’s exact 
test gave the significance value of >3cr. This suggests that, if the companions are separated by 
~280-450 AU, dust emission occurs more frequently in binaries than in single stars. 
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Fig. 7. The figure is the same as Figure 6, but the comparison was performed between closer binary 
systems ( a p < a c ) and single stars with calibration on putative tight companions. Left: Fisher’s p- value vs. 
a c plane. The horizontal-dotted and dashed-dotted lines indicate the 2er and 3er levels, respectively. The 
p-value is minimized at ~450 AU, where the difference in EF between binary systems with a p < a c and 
single stars becomes most significant (2-3er for primary+secondary sources). Right: The EF of binary 
systems and single stars. The horizontal-dotted line indicates the EF s i ng ie (70%). 

However, since the 61 single stars in Ophiuchus have been less studied in high-angular- 
resolution compared to the stars in Taurus, the sample may include undetected close binaries. 
Considering a possible underestimate of EF single , we attempted to calibrate EF single assuming 
the same fraction of close binaries (r/ c b) and the rate of disk-bearing ones (?7 c b;disk) in Ophiuchus 
as in Taurus. Kraus et al. (2012) found that about 30% of their Taurus targets turned out 
to be binaries with 1-40 AU separations, and one third of the binaries showed the signatures 
of disks (rj c b = 0.3, r) c b ; disk = 1/3) . Note that the presence of a disk was inferred through the 
excess emission in 2-8 /irn in their study. If the EF of the single stars in Ophiuchus (hereafter 
EF s;op h) was determined with the same r ] c b and ?fcb;disk, it yielded EF s;op h = 0.659. Thus, the 
number of single stars was estimated as 42 for Ophiuchus, in which 28 were considered to have 
excess emission. Combining the Taurus sample with the calibrated one for Ophiuchus, EF sing i e 
was found to be 70%. We conducted Fisher’s exact test using this sample. The Fisher’s test 
demonstrated that the EF for each primary and secondary of the binaries is still significantly 
higher than that of single stars of 70% if the binary separation is a p ~280-450 AU (~2a; 
Figure 7). The significance of the difference is 2-3cr when estimated for the binary sample 
including both primary and secondary sources (see Figure 7). Therefore, we conclude the 
longer lifetime for binaries in this separation range. 
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6. Discussion 


6.1. EF distribution as a function of disk outer radius 

The binary systems with projected separations, a p , closer than ~450 AU (EF ~ 100%) 
exceeded both wider systems and single stars at the 2-3 a significance level in EF. In binary 
systems, the outer edge of the circumstellar disk is predicted to be truncated and prevented 
from expanding freely by a companion star (e.g., Artymowicz & Lubow 1994). This truncated 
radius, R ou t , strongly depends on the semi-major axis of the binary system, a. On the other 
hand, disk dispersal timescale is probably governed by multiple stellar and disk properties, but 
at least it should be linked to the size of the disk. Therefore, the relationship between EF and 
a p can be interpreted by R ou t . We constructed a model of EF distribution as a function of R 0ut , 
EF mode i(i? out ), to explore whether the characteristic disk size exists to explain the observed high 
EF at the projected separations of 280-450 AU. In the model, we divided the radial regions 
into three sections and assume a constant EF within one section. The radial ranges of the three 
regions and the EFs are treated as variables. To figure out the appropriate ranges of EF mo d e i 
and Rout for this model, we calculated EF mo d e i(a p ) through EF mo d e i(Rout) and compared it with 
the observed one, hereafter denoted as EF obs (a p ). 

The relationship between a and a p is given by a/a p = (1 — e 2 )[l — cos 2 (6 l — </>)sin 2 i] a5 /(l + 
ecos#), where 6, 0, i, and e denote the anomaly, reference angle, inclination, and eccentricity, 
respectively. To derive the a/a p distribution, the eccentricity distribution is required. According 
to recent studies, the eccentricity distribution is relatively uniform (e.g., Raghavan et ah 2010; 
Dupuy & Liu 2011), and high-eccentricity (e> 0.8) binaries are absent. Therefore, we assume a 
uniform eccentricity distribution ranging from e = 0 to 0.7. Using this eccentricity distribution, 
a peak appears at a/a p = 1 in the a/a p distribution. In the previous section, we measured 
EF obs (a p ~ 280-450 AU) = 100% and EF obs (a p > 450 AU) ~ 70%; therefore, we expect that 
EF mo dei(a — 280-450 AU) = 100% and EF mo d e i(a > 450 AU) ~ 70% would give a EF distribution 
consistent with the observations. In addition, EF mode i(a 40 AU) ~ 20-40% is assumed based 
on the report by Cieza et al. (2009). Adopting a truncation radius of R ou t = 0.337a found 
in the systems of mass ratio q = M s /M p = 1 with low-eccentricity (Pichardo et ah 2005), 
the EF can be estimated as a function of R out over three separation regions; R out < 10 AU, 
Rout ~ 100-150 AU, and R ou t > 150 AU. From the above arguments, the EF mo d e i(Rout)’s in 
the three separation regions are varied around ~20, 100, and 70%, respectively. Between the 
three regions, the EF is set by linear interpolation, as shown in Figure 8. The inner and outer 
limits of the radial region of EF mo dei(Rout) — 1 are denoted as R ou t ; i and R out; h, respectively. We 
also define R out;c h and R out;w i as the outer boundary of the EF mode i(Rout) — 0.2 and the inner 
one of the EF mo d e i(R 0 ut) — 0.7 region, respectively. From the relationship between EF mo d e i and 
R out determined from free parameters {R ou t ; i, R ou t ; h) R 0 ut;ch) Rout ; wi} and the EFs in the three 
separation regions, we can derive {a p , R ou t, EF mo d e i(a p )} from randomly selected binary orbital 
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Fig. 8. Schematic of the EF mode i(I? out ) model. This figure shows one example case of three EFs. R a u t;ch 
expresses the outer radius of the innermost EF region, where the EF is set to be similar as observed in 
binaries separated by less than ~40 AU (Cieza et al. 2009). R ou t:i and -R 0 ut ; h represent the inner and 
outer limit of the high EF region, respectively. Introducing these two variables yields a “peak” in the EF 
distribution. R aut;w i represents the inner boundary of the outermost EF region. The three regions are 
linearly connected (dotted lines in the figure). 

parameters The calculation assumes a uniform distribution of semi-major axes. 

This Monte Carlo approach allows us to compare the EF mode i(a p ) and EF obs (a p ) in order to 
determine the parameters related to the EF model (i? out ). 

Figure 9 shows the results of the Monte Carlo computations obtained with R 0 ut; i, R 0 u t ; h, 
-R ou t;ch, and i?out ; wi set to 30-40, 100, 10-20 AU, and 150-250 AU, respectively. The EFs of 
the three characteristic regions are 10% (<R 0 u t ; ch), 100 % (-R 0 ut;m-R out; h), and 70% (>i? 0 ut ; wi)- 
Although the flat EF regions shown in Figure 8 are smoothed by the projection effect, the 
EF remains high around a p ~ 280-450 AU. In other words, the observed high EF region was 
reproduced when introducing this “peak” of EF into the model distribution. The model results 
without peaks are shown in the figure for comparison. For the model without a peak, we set 
EF mode i(>i? out ;i) = 70% so that the EF mode i(a p ) at larger a p approaches the observed EF for 
wider binaries and single stars. We also set EF mode i(i? out;c h) = 10-20% so that the EF of the 
inner 40 AU is consistent with the observed one, 20-40% (Cieza et al. 2009). These results 
disagree with the observational data around a p ~ 280-450 AU (Figure 9). 

The separation-dependent distribution was investigated individually for excess and non¬ 
excess sources, by a Kolmogorov-Smirnov (K-S) test. We iteratively sampled ~10 4 sources from 
the EF mode i(i? out ) at the same a p as the observed binary, and counted the number of excess 
or non-excess objects at each a p . In this way, our model predicts how many sources show 
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Fig. 9. EF moc j e i(a p ) derived from EF mo d e i(-R ou t). Black and gray solid lines denote R 0 u t ;c h = 10 and 
20 AU, respectively. In the peak model, the EF at wide binary separations (a p > 100) depends on R ou t ;w i- 
The selected parameter values are indicated in this figure. If no peak is present in the distribution, the 
EF is almost flat at wide separations. The observational EF data for primary sources are shown for 
comparison. Error bars represent the ler confidence interval estimated from the binomial distribution. As 
reference, the EF of a closely spaced binary (<40 AU, from Cieza et al. 2009) is plotted with dotted error 
bars. The dotted horizontal line represents the EF of a single star (see Section 5.3). 

and do not show excess at a p and their cumulative distribution functions (CDF). According 
to the K-S test, the model predicts the separation-dependent distribution of both primary and 
secondary excess sources, regardless of whether the peak is included in the model. However, 
for sources that do not show excess, the model predictions disagree with observations when 
the peak is excluded. As shown in Figure 10, p >0.32 is obtained by the K-S test where p is 
the probability that the observed binary sources (primary plus secondary) are drawn from the 
same parent sample as the model with an EF peak. On the other hand, p -value is <0.1 for the 
model without a peak, suggesting disagreement with the observations. Therefore, although it 
is not statistically significant, the peak-ed EF mo d e i is more consistent with the observed result 
than the flat EF model . We discuss the interpretation of the EF peak seen in a p ~ 280-450 AU 
and attempt to explain why the EF of this separation range is higher compared to that of the 
single stars in the next section. 

6.2. Mechanism of dust dispersal in binary systems 

The EFs of primary and secondary sources in our sample separated by a p ~ 280-450 AU 
are 100jli 7 % and 9ll® 8 %, respectively. At wider separations, the EF reduces to that of single 
stars (70 ±5%). On the other hand, the EF of close binaries (a p ~ 10-40 AU) has been 
previously determined as ~20-40% (Cieza et al. 2009). In the previous section, we expressed 
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Fig. 10. Cumulative distributions of the projected separation, a p . As demonstrated by the K-S tests, if 
the “peak” is excluded from the EF moc j e i(i? ou t) distribution, discrepancy appears between the predicted 
and the observed separation-dependent EF distribution. Solid and dashed lines show the results with 
and without the peak in the model, respectively. Black and gray colors denote the calculations assuming 
i? out;c h = 10 and 20, respectively. The observed distribution is indicated by the dash-two-dotted line. 

the separation-dependent EF in terms of Rout- The EF significantly decreases, relative to that 
of single stars, when R out < 10 AU, increases when R out ~ 30-100 AU, and matches the single 
star value at R out > 150 AU. This has the implication: the time scale of dust dissipation might 
he prolonged in disks of characteristic radius (~30-100 AU). This disk size is smaller than 
the size estimated for single stars, reported as ~200 AU (e.g., Andrews & Williams 2007) and 
~100-1100 AU (Schaefer et al. 2009 and references therein). Given that the EF is reduced in 
wide binary systems (>450 AU, thus R out > 150 AU), we speculate that extended disks (greater 
than ~30-100 AU) undergo significant mass loss in the outer disk region. 

Various mechanisms of gas and dust removal have been proposed, including viscous 
evolution, planetesimal growth, photoevaporation, disk wind induced by magnetorotational 
instability, and sculpting by planets (e.g., Williams & Cieza 2011; Suzuki et ah 2010). In 
viscous evolution, the timescale of dissipation is given by t oc Rl ut . The R out of binary systems 
depends upon the semi-major axis a, being approximated by ~a/3 (e.g., Artymowicz & Lubow 
1994; Pichardo et al. 2005). Therefore, the peak in the EF distribution is difficult to reproduce 
assuming viscous evolution alone, since more time is required to disperse the dust disk as the 
binary separation (hence the disk size) increases. 

Disk photoevaporation is a well-known mechanism and may successfully explain the 
actual disk dispersal (e.g., Hollcnbach et al. 1994; Clarke et al. 2001; Alexander & Armitage 
2007). Ionizing flux from the central star creates a disk thermal wind and material beyond the 
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critical radius can escape the disk as the gravitational potential weakens. On the other hand, 
viscous accretion continues, leaving an inner hole within the critical radius due to the lack of 
sufficient mass supply from the outer region and the short timescale of viscous process at small 
radii. Armitage et al. (2003) considered the photoevaporation in their study of disk evolution 
in binary systems. In their modeling work, the disk lifetime in a binary depended on its 
separation, and the maximum lifetime was predicted at ~15-100 AU, depending on the model 
which they used. They speculated that binary systems with those intermediate separations have 
less area where the disk mass can be lost by photoevaporative wind and moderately counteract 
the disk dispersal. This earlier prediction cannot be verified by the results presented here 
because we cannot ascertain whether the sources with no IRAC excess really lack outer disks 
without analyzing the longer wavelengths data. However, it is worth pointing out that the EF 
distribution has a peak in the intermediate separation range also in our study. This similarity 
between the observed result and the theoretical prediction implies that photoevaporation has 
a possibility of being one of the effective mechanisms controlling the disk evolution in binary 
systems. 

Disk lifetime is related also to disk mass. If we assume the disk dispersion described 
above, among the binaries of a p > 280 AU, the closer binaries ( a p < 450 AU) should have heavier 
disks than the wider binaries (a p >450 AU). However, such a trend has not been found so far. 
Binaries separated by >280 AU are speculated to have the same degree of millimeter luminosity, 
in other words, the same disk mass as single stars (Harris et al. 2012). Unfortunately, it is still 
difficult to statistically confirm the change in disk mass within a narrow separation range of 
~ 280-500 AU. 

There is another possibility to interpret the dependence of EF on binary separation. 
Systems with wider separations tend to be triples or quadruples. If the 450 AU represent the 
boundary beyond which binaries become sufficiently wide that they start to have additional 
higher-order components, and they harbor close companions, the decreasing of EF can be ac¬ 
counted for by an increasing of existence of multiple systems. In order to clarity this possibility, 
a larger sample is needed for wide binaries for which the presence of close companions were 
thoroughly investigated. 

7. Summary 

By using archival data obtained by Spitzer/ IRAC, infrared excess emission was investi¬ 
gated for 27 binary systems in two nearby star-forming regions, Taurus and Ophiuchus. Our 
sample consisted of binaries with projected separations of 2. ,, 0-17."0 (~280-2400 AU), enabling 
us to study circumstellar disks associated with individual binary components. 

After excluding Class I and objects that cannot be young, spatially resolved photometry 
yielded the following results for 24 binaries or 51 primary and secondary sources. 
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1. The excess frequency (EF) of all primary and secondary stars based on [3.6] — [8.0] color 
is indistinguishable from that of single stars (70 ±5%). This result is consistent with past 
studies in which the disk frequency of wide binary systems (>100 AU) is not significantly 
different from that of single stars. The sources showing excess emission at 8.0 /am also 
exhibits excess at the shorter wavelengths. Conversely, sources without excess at 8.0 /am 
are photospheric at the other IRAC wavelengths. The exceptions are three sources that 
show excess at 8.0 /am but not in one or more shorter wavebands. 

2. In four systems out of 24 (17ig 1 %), excess emission is shown by either the primary or 
secondary component (designated “mixed” systems). In one of these mixed systems, excess 
is found only for the secondary star. No significant preference for either the circumprimary 
or circumsecondary disk to disperse is found. However, the rareness of mixed systems 
suggests that dust dispersal is strongly correlated between the two disks, even at the 
binary separation exceeding >280 AU. 

3. The EFs of closely and widely spaced binary systems are different, consisting of the two 
groups divided by a p ~450 AU. For the entire sample (including both primary and sec¬ 
ondary sources), this difference is significant at the 2-3<r level. The EF of primary and 
secondary sources separated by less than 450 AU was 100l? 7 % and 9llf 8 %, respectively. 
This high EF was compared with that of single stars within a similar range of luminosity 
and spectral type (i.e., similar age and mass). The Fisher’s exact test shows that the 
maximum EF (~100%) is higher than that of single stars with a significance of >2-3 a for 
three kinds of binary sample consisting of primary or secondary or both components. 

Comparing the predicted separation dependence of the EF distribution with the ob¬ 
servational data, the model requires a high EF region (“peak”) around R out ~30-100 AU, as 
confirmed by the K-S test. This disk size is smaller than the typical one estimated for single 
stars, implying that the dissipation of moderately truncated circumstcllar disks is prolonged. 
The EF of widely-spaced binary systems, presumably with larger circumstcllar disks, is indis¬ 
tinguishable from that of single stars. Our result is consistent with the prediction that extended 
disks may suffer from significant mass loss in the outer regions by photoevaporation (Armitage 
et al. 2003), suggesting that photoevaporation contributes to the evolution of circumstellar 
disks in binary systems. 

We are grateful to the referee for valuable comments to improve the paper. This 
work is based on archival data obtained by the Spitzer Space Telescope, which is operated 
by the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the 
National Aeronautics and Space Administration (NASA). Support for this work was provided 
by NASA. This publication utilizes data products from the Two Micron All Sky Survey, which 
is a joint project of the University of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology and is funded by the NASA and the National Science 
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Appendix 1. Notes on individual binaries 

A.1.1. 2MASS J16263682-2415518 B 

This binary system is located in Ophiuchus. The secondary star begins to exhibit signif¬ 
icant excess (0.48 ±0.14 mag) at 8.0 /irn. The binary separation is very wide at 9"08 (Duchene 
et al. 2007), and the spectral type of the secondary star is estimated as M5 (Wilking et al. 
2005). 

The SED fitting was performed to the flux density in R , I, and J bands obtained 
from USNO-B1 and 2MASS Catalog. Significant excess emission was confirmed at 8.0 /ini, 
but was not identified at 5.8 /irn. However, the possible excess at 5.8 /ini might have been 
underestimated because the photometric values at 3.6 and 4.5 /irn are smaller than the fitted 
photosphere. The measured color of [4.5] — [5.8] (0.49 ±0.2 mag) is in fact significantly higher 
than the photospheric color for an M5-type star (0.04 mag; Luhman et al. 2010), thus excess 
emission likely begins at 5.8 /im. The discrepancy between the photometry and the fitted 
photospheric values at 3.6 and 4.5 /im could arise from the poor SED fitting due to the faintness 
of the object in the R band (19.2 mag). Therefore, according to its color, we consider that this 
object has an onset of the excess at ~5.8-8.0 /irn. 

The source was detected at 24 /im in the c2d survey, and the photometric value of 
5.71 ±0.19 mag was provided in the catalog. The [8.0] — [24] color is then 3.46 ±0.24 mag, which 
is significantly higher than the photospheric one for an M5-type star (0.23 mag; Luhman et al. 
2010). Thus, we again confirmed the existence of circumsecondary dust, and the occurrence of 
dust dispersion in the inner region of the disk. 

A. 1.2. VllOTauA 

The binary separation of V710 Tau is 370 (White & Ghez 2001). The spectral types of 
the primary and secondary stars are estimated as M0.5 and M2, respectively (White & Ghez 
2001). In V710 Tau, excess emission is found only for the primary source. The subtraction of 
the photospheric model derived from the SED fitting resulted in significant excess at 5.8 and 
8.0 /irn, with values of 0.37 ±0.14 and 0.94 ±0.14 mag, respectively. On the other hand, no 
excess was observed at 3.6 and 4.5 /im. The [4.5] — [5.8] color of 0.67 ±0.2 mag is significantly 
higher than the photospheric color of 0.2 mag, also suggesting that excess emission begins from 
5.8 /irn. In other wavebands, the primary star shows a small excess in L band ( K — L = 0.44; 
White & Ghez 2001) and a substantial one in N (K — N = 3.32; McCabe et al. 2006). The 
significant excess starting from 5.8 /mi is not inconsistent with these previous studies. 

The flux density at 24 /im was obtained as 236 mJy for the whole system (primary plus 
secondary) (Wahhaj et al. 2010). The magnitudes measured at 8.0 and 24.0 /un for the system 
are 6.30 ±0.20 mag and 3.71 ±0.11 mag, respectively, then the color of [8.0] — [24] is 2.59 ± 
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0.23 mag. This is significantly larger than an averaged-photospheric color of an M-type star of 
0.22 mag (Luhman et al. 2010). Because the emission at 24 /irn cannot be spatially resolved, it 
is difficult to conclude whether the origin of the emission is circumprimary, circumsecondary, 
both of them, or circumbinary disks. However, considering the separation of ~ 280 AU, the 
emission at 24 //m is possibly not from the inner edge of a circumbinary disk. 

A.1.3. JH 223 A 

The binary JH 223 belongs to the Taurus region. The primary star of this system shows 
the significant excess at 4.5 //in. The binary separation is small (2707; Kraus et al. 2011), and 
the spectral types of its primary and secondary are estimated as M2 and M6.5, respectively 
(Kraus et al. 2007). 

The excess emission at 3.6, 4.5, 5.8, and 8.0 //m were —0.04 ±0.03, 0.33 ±0.04, 0.61 ± 
0.04, and 1.03 ±0.04 mag, respectively, estimated by subtracting the photospheric magnitudes 
obtained from the SED fitting. This indicates that the onset of the excess is at around 4.5 fi m. 
The resultant x 2 /d.o.f. values for the PSF fitting photometry were 1.0, 0.4, 0.7 and 1.0 at the 
four wavelengths, respectively, implying that the fitting was reasonably performed. The Ha 
emission have identified this source as a Weak-lined T Tauri Star (Neuhaeuser et al. 1995 and 
references therein), suggesting that this source lacks an optically-thick, inner disk with active 
accretion, but may possess an optically-thick, exterior disk. 

The photometry at 24 //m was previously obtained as 5.13 ±0.04 mag for the primary 
and secondary sources (Rebull et al. 2010). The measured magnitudes at 8.0 //m of the system is 
7.72±0.20 mag and the [8.0] — [24] color is calculated to be 2.59±0.20 mag, which is significantly 
higher than an averaged-photospheric value of an M-type star of 0.22 mag (Luhman et al. 2010). 
As in the case of V710 Tau, although it is hard to conclude the origin of the excess emission, 
the emission at 24 /xm may not be caused by the circumbinary disk. 
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